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I  INTRODUCTION  AND  BACKGROUND 


The  U.S.  Naval  Observatory  in  Washington,  D.C.,  supplies  an  atlas 
of  the  absolute  positions  of  celestial  bodies  to  within  .01  seconds  of 
angle  for  astronomical  research  and  navigation.  (These  are  fundamental 
measurements;  the  relative  positions  of  other  heavenly  bodies  are  not 
used.)  A  major  technique  used  for  this  work  is  Transcircle  Observation. 
Employing  a  combination  of  the  Transit  Circle  Telescopic  instrument  with 
integral  north-south  pivoting  and  the  earth's  east-west  rotation  during 
observations  allows  full-sky  Leverage. 


A.  Problem  Description 

Within  the  last  few  years, improvements  in  instrumentation  precision 
(both  mechanical  and  optical)  have  led  to  the  realization  that  the  re¬ 
moval  of  systematic  errors  due  to  meteorological  causes  (often  with  large 
spatial  and  temporal  variations)  is  a  major  consideration  in  stellar 
position  accuracy  and  in  the  repeatability  of  the  arrangement  of  posi¬ 
tions  taken  during  different  nightly  observations.  The  effects  of  the 
atmosphere  tend  to  be  site-specif ic ,  and  perhaps, include  seasonal  varia¬ 
tions.  Thus,  attempts  to  handle  the  variation  in  measured  positions  with 
multiple  sightings  from  the  same  observatory  location  and  then  "averaging" 
the  error  to  reduce  its  magnitude  will  not  improve  the  systematic  in¬ 
fluence  on  the  error  structure  of  the  estimates. 

Although  each  major  atmospheric  constituent  can  affect  viewing  the 
stars,  two  major  areas  have  been  singled  out  as  having  the  highest 
priority  (see  the  Appendix).  The  first  is  water  vapor  content  along 
the  viewing  path.  Water  vapor  changes  the  gas  constant  and  refractivity 
values  (from  those  measured  and  used  when  assuming  dry  air).  Presently, 
the  ground-level  value  of  the  partial  pressure  of  water  and  a  theoretical 
model  for  the  decrease  in  water  vapor  density  with  altitude  are  used  to 
compensate  for  the  optical  bending  applied  to  the  viewing  data.  An 
estimate  of  the  actual  water  vapor  distribution  with  an  accuracy  of 
approximately  20  percent  would  routinely  correct  position  errors  of 
.01  seconds  of  angle.  For  those  worst-case  situations  where  the  existing 
structure  grossly  deviates  from  the  linear  model,  errors  approaching  .1 
seconds  of  angle  could  be  corrected. 

The  second  problem  involves  isopycnic  tilts  (layers  of  equal  atmo¬ 
spheric  density).  An  atmospheric  model  that  entails  the  existence  of 
concentric  spheres  of  equal  density  layers  is  incorrect  and  leads  to 
errors  in  position  when  applied  to  refraction  models,  especially  when 
viewing  at  slant  angles  instead  of  directly  overhead.  The  tilt  (from 
exactly  concentric)  is  typically  a  few  degrees.  Assuming  that  isotherms 
(lines  of  constant  temperature)  follow  isopycnics,  planes  of  constant 
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temperature  can  be  measured  to  specify  the  degree  of  tilt  along  the  view¬ 
ing  path.  Thus,  knowing  the  temperature  structure  with  an  accuracy  of 
.5°K  would  remove  the  systematic  errors  caused  by  this  effect,  which 
range  from  .01  up  to  .1  seconds  of  angle  at  extreme  zenith  viewing 
angles.  Because  of  the  precision  needed  in  atmospheric  modeling  used 
by  the  correction  algorithms,  absolute  measurements  of  the  temperature 
and  water  vapor  content  are  far  more  useful  than  relative  measurements. 


B.  Lidar  as  a  Solution 


Lidar  (Light  Detection  and  Ranging)  provides  a  remote-sensing  solu¬ 
tion  to  the  measurement  problems  of  the  Observatory.  Since  the  measure¬ 
ments  are  a  set  of  planes  at  different  heights  (to  a  maximum  of  5  km) , 
three-dimensional  coverage  with  a  conical  scan  lidar  is  necessary.  A 
plane  can  be  defined  by  measurements  taken  at  different  positions  at 
the  same  altitude.  The  approach  is  similar  to  conventional  radar,  except 
the  interactions  occur  at  optical  frequencies.  With  a  pulsed  laser 
transmitter  emitting  a  narrow,  well-collimated  beam,  and  a  receiver 
system  (tuned  to  respond  to  a  specific  set  of  optical  wavelengths)  posi¬ 
tioned  at  a  common  location,  a  telescopic  optical  receiver  can  perform 
ranging  by  calculating  the  time  for  a  pulse  of  radiation  to  travel  the 
distance  to  a  region  of  interest  plus  the  time  for  the  backscattered 
pulse  to  return.  By  sampling  the  time-varying  received  optical  intensity 
after  a  specific  delay,  study  of  a  restricted  location  along  the  path  can 
be  accomplished. 

The  lidar  system  should  be  a  mobile  unit  so  measurements  could  be 
obtained  at  other  observatory  locations,  and  for  greatly  increasing  its 
value  as  a  general  development  tool. 

For  most  of  the  compensation  programs  used  at  the  Naval  Observatory, 
a  complete  atmospheric  profile  derived  once  an  hour  is  sufficient.  The 
profiles  should  be  measured  and  collected  by  an  integral  data  collection 
computer  system  with  storage  capability  (magnetic  tape  is  the  most  use¬ 
ful)  for  post-measurement  analysis. 

A  final  requirement  and  consideration  is  that  operation  during  the 
daytime  should  not  be  excluded.  Although  not  a  major  requirement, 
important  position  measurements  can  be  taken  during  daylight  hours. 
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II  APPROACH 


Two  laser  systems,  each  using  different  physical  principles,  can 
supply  the  atmospheric  data  needed  by  the  Naval  Observatory.  The  first 
is  Raman  Lidar;  the  second,  Differential  Absorption  Lidar,  or  DIAL 
(Hinkley,  1976).  This  section  will  briefly  describe  each  system. 


A.  Raman  Lidar 


This  system  transmits  a  single  wavelength,  X>p»  anc*  measures  the  back- 
scattered  energy  received  at  up  to  three  other  wavelengths,  Xjq,  XR2>  and 
Xj^3,  all  different  and  uniquely  Raman-shifted  (from  the  transmitted  wave¬ 
length)  by  the  molecular  species  of  interest  (Cooney,  1970).  The  amount 
of  backscatter  received  is  proportional  to  the  Raman  backscatter  coef¬ 
ficient  for  the  molecule  being  excited  at  the  transmitted  wavelength. 
Figure  1  shows  the  change  in  backscatter  from  nitrogen  as  the  ambient 
temperature  varies.  Temperature  estimates  are  calculated  by  comparing 
the  differential  backscatter  measurements  at  Xjq  and  XR2  for  a  change  in 
temperature,  AT  =  T2  -  T]_.  Only  the  relative  temperature  differential 
between  two  heights  can  be  directly  obtained.  By  measuring  the  relative 
backscatter  strength  from  water  vapor  at  XR3,  the  relative  number  density 
of  water  vapor-to-nitrogen  can  be  evaluated. 


FIGURE  1  TEMPERATURE  DEPENDENCE  OF  RAMAN  BACKSCATTER 
FROM  NITROGEN  AS  A  FUNCTION  OF  WAVELENGTH 
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B. 


DIAL  Lidar 


A  DIAL  lidar  system  transmits  at  least  two  wavelengths,  AT^  and  A^j 
tuned  to  correspond  to  an  absorption  line  of  the  species  of  interest, 
the  "on"  wavelength,  and  to  a  nonabsorbing  region,  the  "off"  wavelength. 
(Measurement  techniques  for  temperature  include  a  third  "on"  wavelength.) 
The  receiving  system  measures  the  elastic  backscattered  intensity  at  each 
of  the  transmitted  wavelengths,  i.e.  Ar^  =  ^T1  an<*  ^R2  =  ^T2*  Backscatter 
is  due  to  gaseous  molecular  components  (Rayleigh)  and  to  aerosol  content 
(Mie)  (Browell  et  al.,  1981;  Shotland,  1974). 

The  ratio  in  the  magnitudes  of  the  returned  signals  at  the  "on" 
and  "off"  wavelengths  can  be  inverted  to  yield  a  concentration  estimate. 

A  temperature  estimate  is  obtained  by  using  another  "on"  wavelength, 

A-J3.  By  analyzing  the  variation  of  the  change  in  absorption  strength 
with  respect  to  temperature  at  the  two  "on"  wavelengths,  the  absolute 
temperature  can  be  measured.  The  selection  of  the  "on"  absorption  region 
is  a  tradeoff  between  choosing  a  line  of  sufficient  line  strength  such 
that  the  differential  measurements  can  be  obtained  with  sufficient 
signal-to-noise  qualities,  and  yet  not  have  the  absorption  of  the  trans¬ 
mitted  laser  energy  be  so  great  as  to  prevent  enough  backscattered  signal 
to  be  received  from  the  maximum  desired  range.  Section  III.C  details 
this  technique;  Section  III.D  describes  transmission  characteristics  at 
735  nm,  an  optimum  wavelength  region  for  laser  transmission. 


C.  Water  Vapor  Profiling  Systems 

Both  the  Raman  and  DIAL  lidar  techniques  have  successfully  demon¬ 
strated  vertical  water  vapor  density  measurements.  Significant  experi¬ 
ments  that  have  supplied  vertical  profiles  are  reviewed  in  Tables  1 
and  2.  Those  systems  that  obtained  a  column- Integrated  measurement  or 
an  average  value  for  a  horizontal  path  were  excluded. 

Table  1  lists  Raman  lidar  systems  that  have  performed  vertical 
water  vapor  profiling.  All  these  experiments  measure  the  ratio  of  the 
number  desity  of  water  vapor  to  that  of  nitrogen.  Thus,  only  a  water 
vapor  mixing  ratio  is  estimated.  (This  can  be  translated  into  an  absolute 
water  vapor  density  profile  by  assuming  a  model  density  profile  for 
nitrogen  with  accompanying  error.) 

Table  2  lists  DIAL  systems  for  vertical  water  vapor  profiling  Less 
powerful  laser  transmitters  are  needed  than  those  used  in  Raman-based 
systems  because  the  Rayleigh-Mie  backscatter  cross  sections  are  many 
orders  or  magnitude  larger  than  Raman  cross  sections.  (Note  that  the 
DIAL  systems  have  been  operated  during  the  daytime,  whereas  Raman  sys¬ 
tems  have  been  limited  to  night  operation  only.) 
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Table  1 


RAMAN  LIDAR  SYSTEMS  FOR  VERTICAL  WATER  VAPOR  PROFILING 

Cooney  (1970) 

Transmitter:  Doubled  ruby  at  347.15  nm. 

Receiver:  397.5  nm  from  H^O  and  377.7  nm  from  N£. 

Approach:  Water  vapor  prof^xes  obtained  by  normalizing  H2O  Raman  re¬ 
turn  with  N2  return  and  assuming  an  N2  atmospheric  vertical 
model  content. 

Accuracy:  13  %. 

Range:  300  m  -  2600  m  with  100-m  intervals. 

Nighttime  measurement. 

Melfi  (1972) 

Transmitter:  Doubled  ruby  at  347.15  nm. 

Receiver:  397.5  nm  from  H^O  and  377.7  nm  from  Nj. 

Approach:  Water  vapor  profile  from  H2O  return  normalized  by  N2 
return. 

Range:  200  m  -  2500  m  with  50-rn  intervals. 

Nighttime  measurement. 

Pourny,  Renault,  and  Orszag  (1979) 

Transmitter:  Doubled  ruby  at  347.15  nm. 

Receiver:  397.6  nm  from  H2O  and  377.7  nm  from  N^. 

Approach:  Water  vapor  mixing  ratio  from  H^O  return  normalized  by  N2 
return. 

Accuracy:  15% 

Range:  300  m  -  1800  m  with  30-m  intervals. 

Nighttime  measurement. 

Renault,  Pourny,  and  Capitini  (1980) 

Transmitter:  Quadrupled  YAG  at  366  nm. 

Receiver:  277.5  nm  from  O2,  283.6  nm  from  N2,  and  294.6  nm  from  H2O. 

Approach:  Water  vapor  mixing  ratio  from  the  H2O  return  and  is  normal¬ 
ized  by  the  N2  return.  Estimation  of  ozone  concentration 
is  used  to  compensate  for  heavy  ozone  attenuation  of  the 
signals  in  thi3  wavelength  region. 

Accuracy:  10%  at  500  m  with  30-m  intervals. 

Range:  150  m  -  950  m  with  30-m  intervals. 

Daytime  operation. 
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Table  2 


DIAL  LIDAR  SYSTEMS  FOR  VERTICAL  WATER  VAPOR  PROFILING 

Browell,  Wilkerson,  and  McElrath  (1979) 

Transmitter:  Ruby  at  694.3  nm  and  dye  (pumped  by  ruby)  at  724.3  nm. 

Receiver:  724.3  nm  (water  vapor  absorbing  wavelength)  and  694.3  nm. 

Approach:  Difference  in  returned  signals  is  analyzed  for  H2O  density 

profile  (mol/ cm3).  Equipment  includes  an  absorption  cell 
that  measures  the  cross  section  for  H^O  at  the  dye  laser 
output  for  each  pulse. 

Accuracy:  9%. 

Range:  100  m  -  3000  m  with  100-  to  180-m  intervals. 

Nighttime  operation. 

Cahen,  Pelon,  Flamant,  Lefere,  Chanin,  and  Megie  (1980) 
and  Cahen  and  Lesne  (1981) 

Transmitter:  Nd:YAG  at  532  nm  and  dye  (pumped  by  Nd:YAG)  at  723  nm. 

Receiver:  723  nm  (absorption  by  water  vapor)  and  532  nm. 

Approach:  Different  returns  to  obtain  water  vapor  profiles.  Wave¬ 
length  and  line  width  of  transmitters  are  servo-controlled. 

Accuracy:  10%. 

Range:  8000  m  with  30-m  intervals. 

Daytime  operation. 

Werner  and  Herrmann  (1981) 

Transmitter:  Ruby  at  694.28  nm  and  temperature-varied  ruby  from 
694.20  nm  to  694.28  nm. 

Receiver:  694.215  rm  (absorbed  by  water  vapor)  and  694.237  nm. 
Approach:  Difference  returns  to  obtain  water  vapor  density  profiles. 

Accuracy:  ^  rabar  in  absolute  humidity. 

Range:  450  m  -  1500  m  with  100-m  intervals. 

Daytime  operation. 
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D.  Temperature-Profiling  Systems 


Table  3  contains  a  list  of  the  Raman  lidar  systems  for  vertical 
temperature  profiling.  Common  to  these  experiments  is  the  tchnioue  of 
measuring  the  temperature  dependence  in  the  magnitude  of  the  Raman  spec¬ 
trum  of  oxygen  and  nitrogen.  The  lines  in  the  Raman  spectrum  :  r  these 
molecules  have  differential  changes  in  magnitude  with  respect  to  tempera¬ 
ture.  Thus,  by  comparing  the  differential  change  in  two  bands,  a  tempera¬ 
ture  value  can  be  obtained. 

There  is  a  lack  of  DIAL  vertical  temperature  profiling  svstems; 
however,  work  is  presently  being  performed  with  techniques  that  use  water 
vapor  lines.  One  approach  is  to  use  "Mason's  Method"  with  the  absorption 
spectra  of  water  vapor  in  the  region  of  735  run  (see  Section  III. C).  This 
wavelength  is  close  to  those  used  for  successful  DIAL  water  vapor  mea¬ 
surements  (724  nm)  and  other  temperature  techniques  using  the  absorption 
spectra  of  oxygen  (760  nm)  (Kalshoven  et  al.,  1981).  Part  of  the  interest 
in  the  area  has  been  generated  by  the  introduction  of  a  tunable  solid- 
state  alexandrite  laser  that  has  direct  application  for  study  in  this 
wavelength  region  (see  Section  III.D). 
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Table  3 


RAMAN  LIDAR  SYSTEMS  FOR  VERTICAL  TEMPERATURE  PROFILING 


Cooney  and  Pina  (1976) 

Transmitter:  Ruby  at  694.3  nm. 

Receiver:  689.0  nm  and  691.2  nm. 

Approach:  Measure  the  difference  in  the  magnitudes  of  the  Raman  spec¬ 
trum  of  O2  and  N2  in  two  bands  which  vary  with  temperature. 

Accuracy:  5°C. 

Range:  600  m  -  1100  m  with  300-m  intervals. 

Nighttime  operation. 

Gill,  Geller,  Farina,  and  Cooney  (1979) 

Transmitter:  Ruby  at  694.3  nm. 

Receiver:  691.5  nm  and  689.0  nm. 

Approach:  Raman  spectrum  of  O2  and  N2* 

Accuracy:  .85°C. 

Range:  1300  m  -  2300  m  with  75-m  intervals. 
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Ill  ANALYSIS 


Our  evaluation  of  the  capabilities  of  each  lidar  technique  (DIAL 
versus  Raman)  as  well  as  evaluations  by  other  researchers  (Wright  et  al., 
1974;  Hake  et  al. ,  1977)  indicate  that  a  DIAL  system  should  be  directly 
applied  to  the  measurements  needed  by  the  Observatory.  The  DIAL  tech¬ 
nique  has  demonstrated  sensitivity  for  water  vapor  content  profiling  and 
can  be  designed  to  provide  sufficient  atmospheric  temperature  structure 
information  during  daytime  operation. 

The  study  of  the  DIAL  technique  involves  the  use  of  absorption 
spectra  for  constituent  and  trace  gases  in  the  atmosphere  and  numerical 
models  for  the  transmission  of  optical  beams  through  the  atmosphere. 
Research  and  the  collection  of  experimental  and  theoretical  data  by  the 
Air  Force  Geophysics  Laboratory  (AFGL)  have  provided  widely  recognized 
code  (computer  programs)  on  optical  propagation  (McClatchey  et  al.,  1970). 
To  use  this  compendium  for  numerical  analysis  in  optical  propagation, 
code  that  describes  the  absorption  of  laser  lines  in  the  atmosphere  had 
to  be  converted  from  the  computer  language  Fortran  IV  for  Control  Data 
Corporation  computers  to  Fortran  77,  compatible  code  for  use  on  the 
Digital  Equipment  Corporation  VAX  11/782  installed  at  SRI  International. 


A.  Absorption  Spectra  Compilation 

The  absorptic  spectra  is  supplied  on  two  magnetic  tapes  available 
from  AFGL.  One  is  called  the  Major  Constituents  Compilation  (Rothman, 
1980;  McClatchey,  1973);  the  other  is  the  Trace  Gas  Compilation  (Rothman 
et  al. ,  1981).  Files  on  the  major  constituents'  tapes  contain  the 
molecular  transitions  of  the  seven  major  gases  in  the  atmosphere: 

•  Water  vapor 

•  Carbon  dioxide 

•  Nitrous  oxide 

•  Carbon  monoxide 

•  Methane 

•  Ov^rfron 

J 

•  Ozone. 

The  total  number  of  absorptions  lines  is  over  159,000,  with  the  spectra 
ranging  from  .3  to  17,880  cm~l  (560  nm) . 

The  absorption  lines  of  13  trace  gases  are  collected  in  a  separate 
atlas;  these  include: 
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•  Nitric  oxide 

•  Sulfur  dioxide 

•  Nitrogen  dioxide 

•  Ammonia 

•  Nitric  acid 

•  Hydroxyl  radical 

•  Hydrogen  halides  (HF,  HC1,  HBr ,  HI) 

•  Chlorine  monoxide 

•  Carbonyl  sulfide 

•  Formaldehyde. 

In  the  study  of  the  DIAL  technique  for  water  vapor  and  temperature 
measurement,  particular  attention  was  paid  to  the  700-  to  800-nm  region 
in  the  infrared  spectrum  where  numerous  water  vapor  transitions  exist. 
This  is  also  a  spectral  region  relatively  free  of  absorption  lines  from 
other  molecules. 


B.  Computer  Codes 

The  computer  program  FSCDATM  represents  the  earth's  atmosphere  as 
layers  in  thermal  equilibrium.  A  number  of  atmospheric  model  profiles 
can  be  chosen  as  input: 

•  Tropical  atmosphere 

•  Midlatitude  summer  or  winter 

•  Subarctic  summer  or  winter 

•  U.S.  1962  standard. 

The  user  of  the  code  specifies  the  wavelength  region  of  interest 
and  the  geometry  of  the  transmission  path,  the  altitude  of  transmitter 
and  receiver,  and  the  slant  angle  between  them.  The  code  then  outputs 
a  layered  model  of  the  atmosphere  giving  the  total  density,  the  average 
temperature,  and  the  molecular  column  densities  for  each  layer.  The 
layer  thickness  is  varied  according  to  the  user-specified  maximum  tempera¬ 
ture  difference  between  layers.  (The  layer  temperature  differential  di¬ 
rectly  affects  the  computation  speed  and  required  accuracy  of  the  absorp¬ 
tion  calculations.) 

The  derived  model  for  the  atmosphere  and  the  path  geometry  are  then 
input  to  the  code  FASCODl-Fast  Atmospheric  Signature  Code  (Smith  et  al., 
1978).  This  code  performs  a  line-by-line  absorption  calculation  using 
a  Voight  lineshape,  a  combination  of  Lorentian  (pressure)  broadening 
that  best  describes  transitions  in  the  lower  atmosphere;  and  Doppler 
(velocity)  broadening,  with  .01  cm“l  accuracy  (Clough  and  Kneizys,  1979). 
All  the  lines  from  the  seven  atmospheric  molecules  and  the  13  trace 
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bases  within  the  wavelength  region  of  interest  are  applied  in  the  laver- 
by-laver  calculation. 

The  total  path  transmission  is  obtained  by  first  obtaining  the  laver 
parameters  from  FSCDATM.  FASC0D1  then  calculates  the  absorption  for 
each  layer,  merges  the  absorption  results  for  all  layers,  and  convolves 
the  derived  atmospheric  absorption  with  a  modeled  instrument  transmission 
function.  (Specific  laser  lines  can  be  calculated  by  assuming  a  zero 
bandwidth  for  the  laser  source.) 

In  the  next  section  we  describe  a  DIAL  technique  where  the  transmis¬ 
sion  characteristics  of  water  vapor  absorption  lines  can  measure  range- 
resolved  values  of  water  vapor  density  and  temperature  simultaneously. 
Section  III.D  contains  results  of  transmission  simulations  that  verify 
this  technique. 


C .  A  DIAL  Technique  for  Simultaneous  Water  Vapor 

and  Temperature  Measurements 

This  approach  uses  three  wavelengths:  XQ,  X]_,  and  X 2 ,  where  \t,  is 
located  in  a  ncnabsorbing  part  of  the  water  vapor  spectrum,  and  X^  and  Xi 
are  located  at  the  center  of  two  absorption  lines  of  water  vapor  with 
similar  absorption  coefficients  but  different  line  strength  sensitivities 
with  respect  to  temperature  (Endemann  and  Byer,  1981;  Endemann  and  Bver, 
1980;  Lebow  et  al. ,  1982;  Rosenberg  and  Hogan,  1981).  (See  Section  III.D 
for  examples  of  three  wavelengths  near  734  nm.) 


A  DIAL  operation  is  performed  with  the  pairs  (XQ,  X^)  and  p0,  X2). 
The  difference  in  the  optical  depth,  x,  for  one  pair  of  lines  over  a 
range  cell  AR  at  range  R  can  be  written  as: 


x 


1 


‘P1(R)  ' 

P1  (R+AR) 


r  p  (r)  1 

o _ 

P  (R+AR) 


=  2N(R)o1AR  ,  (Hinckley,  1976)  (1) 


where  N(R)  is  the  average  water  vapor  number  density  over  the  cell  AR, 
and  o±  is  the  differential  absorption  coefficient  for  wavelengths  XD 
and  X^.  Atmospheric  backscatter  and  attenuation  caused  by  scattering 
(gaseous  and  particulate)  and  absorption  by  other  species  besides  water 
vapor  are  assumed  to  have  low  spectral  variability  at  the  two  wavelengths. 
A  similar  equation  can  be  written  for  the  second  pair  of  wavelengths. 
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In 


P2(R) 


P2(R+AR) 


-  In 


P  (R) 
o 


P  (R+AR) 
o 


=  2N(R)a2AR 


(2) 


The  differential  absorption  coefficient  is  a  function  of  the  pres¬ 
sure  and  temperature  at  the  range  R.  It  is  defined  as 


o(T,P)  = 


S  (T) 


TT  a(T,P) 


(3) 


where  S(T)  is  the  line  intensity  and  a(T,P)  is  the  line  halfwidth.  In 
general,  the  temperature  dependence  for  the  line  intensity  from  a  nominal 
temperature,  TQ  =  296  K,  is  written  as 


S(T)  - 


S  (T  ) 
o 


Q  (T  ) 
v  o 

Q  (T) 

V 


W 

Q  <T> 


exp 


_k  (t0  -  T  j 


(McClathy ,  1973) 


(4) 


where  Qv  is  the  vibrational  partion  function;  for  water  vapor  Qv  =  1 . 

Qr  is  the  rotational  partition  function;  for  water  vapor  Qr(T)  = 

(T/T0)3/2  Qr (t0) .  E  is  the  energy  of  the  lower  state  for  the  transition, 
and  k  is  Boltzman's  constant.  Then, 


S  (T) 


S (T  )[T(R)/T  ]“3/2 
o  o 


Jj) 

ToJ) 


(5) 


For  the  lower  atmosphere,  the  halfwidth  temperature  and  pressure  de¬ 
pendence  is 


a(T,P)  =  a(T  ,P  )  [T  /T(R)]n  (P/P  )  ,  (6) 
o  o  o  o 

where  PQ  is  a  standard  pressure  for  measurements,  1  atm  =  1016  mbar,  and 
n  is  an  exponential  temperature  factor  that  is  modeled  as  having  a  value 
of  .5. 


The  DIAL  equation,  Eq.  (1) , 
the  temperature  dependence: 


can  then  be  rewritten  explicity  to  show 


P1(R)/P  CR+AR) 

P  (R)/P  (R+AR) 
o  o 


2N(R)AR  • 


s.  (T  )  .  f  T 

1  p 

(e,  r 

loll 

o  0 

)  1 

a.(T  )  tt  |T( 
1  o 

R)J  P  6XP 

) k 

T (R)J  J 


(7) 
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Let  P(R+AR)/P(R)  be  defined  as  the  transmittance,  T (AR) ,  through  the 
cell  AR.  Then, 


’t  (ar)1 
T1  (AR)J  =  2  N 


(R)  AR 


with  a  second  similar  equation  for  the  other  pair  of  wavelengths.  An 
expression  for  the  average  temperature  in  the  range  cell  can  be  written 
by  equating  the  derived  number  density  N(R)  for  each  pair  of  wavelengths: 

ln[To/T1]  S1(TQ)/a1(To)  exp  [E^kTj  exp  [-E^kT] 
ln[T,/T0]  =  S2(To)/a2(To)  exp  [E0/kT_]  exp  [ -E0/kT ]  *  (9) 


Solving  for  T(R), 


<  p  _  F 
2  1 


In  [T  (AR) /T.  (AR) ] 
_ o _ i _ 

In  [T  (AR)  /T  (AR)  ] 
o  2 


-  ln  c. 


and  C2  are  constants  related  to  the  absorption  lines  at  X]  and  A 2  » 
and  can  be  evaluated  using  line  compilation  tables  (e.g. ,  McClatchey,  1973) 
or  laboratory  measurements  (Wilkerson  et  al.,  1979),  i.e., 


Ci  =  [Si(To)/ai(To)]  6XP  [Ei/kT0)] 


Thus,  by  measuring;  the  transmittance  through  a  range  cell  AR  at  three 
wavelengths,  TQ(AR),  T^(AR),  and  T£(AR),  the  average  temperature  for  the 
cell  T(R)  can  be  calculated. 

After  obtaining  the  temperature  for  the  range  cell,  the  water  vapor 
number  density  can  be  calculated  from  the  transmittance  difference  of 
either  of  the  line  pairs.  For  (X0,A^), 


N  (R)  «  ln 


To(AR)  TT  Ctl(T0)  T  (R)  P  (R) 

T, (AR)  2AR  S,  (T  )  T 
1  1  o  o 


exp  hr 


1  (\  ■  t<r>) 


•  (12) 


Since  the  linewidth  is  a  function  of  pressure,  the  accuracy  of  N(R)  is 
dependent  on  an  assumed  value  for  P(R).  However,  by  using  the  ideal  gas 
law,  PV  -  rtRT,  an  accurate  representation  for  the  water  vapor  mixing 
ratio  independent  of  P(R)  can  be  written: 


■To(AR)l 


T1(AR)  2AR  S 


aifV  Tiw  «_  ffi  / 1 _ i_\ 

.  (T  )  T  P  P  k  T  T  (R) 

1  o  o  o  l\°  /. 
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where  ft  is  the  ideal  gas  constant,  and  q(R)  is  the  average  water  vapor 
mixing  ratio  in  the  range  cell  AR. 


D.  Transmission  Characteristics 


Using  the  AFGL  line  compilation  data  and  the  computer  codes  FSCDATM 
and  FASC0D1 ,  numerical  calculations  of  the  transmission  for  radiation 
in  the  region  of  735  nm  were  performed.  This  work  assumes  atmospheric 
models  to  simulate  the  principles  used  in  the  DIAL  technique.  Figure  2 
shows  the  one-way  transmission  spectrum  from  13,600  to  13,630  cm~l  through 
a  vertical  300-m  column  of  1962  U.S.  Standard  Atmosphere  (at  284  K) .  The 
absorption  in  this  region  is  due  solely  to  water  vapor  lines.  The  DIAL 
temperature  measurement  technique  described  in  the  last  section  utilizes 
two  of  the  lines  in  the  region,  13,602  cm~l  and  13,628  cm--*-. 

Figures  3  and  4  demonstrate  the  temperature  dependence  in  the  water 
vapor  lines  in  this  particular  spectrum. 

Individually,  Figure  5  shows  the  change  in  the  transmission  profile 
with  respect  to  temperature  for  the  line  at  13,602.63  cm~l;  Figure  6, 
for  the  line  at  13,628.88  cm"-'-.  Note  that  the  temperature  dependence  of 
these  two  lines  is  inversely  related.  Although  the  absorption  strength 
of  the  two  lines  is  comparable  a  the  standard  measurement  temperature 
of  296  K,  dissimilar  ground-state  energies  for  each  of  the  lines  (920.169 
cm“l  for  13,602.631  and  23.794  cm~l  for  13,628.866)  implies  a  large  dif¬ 
ferential  change  in  absorption  with  an  increase  or  decrease  in  ambient 
temperature. 

Figure  7  represents  the  change  in  transmission  at  line  center  for 
three  wavelengths  when  the  temperature  in  a  range  cell  of  300  m  varies 
from  254  to  314  K.  As  an  implementation  of  the  DIAL  technique  described 
in  Section  III.C,  AQ  (the  reference  wavelength)  is  chosen  to  be  734.4  nm; 
A^  and  A2  would  then  correspond  to  the  two  absorbing  lines  at  733.9  and 
735.1  nm. 


E.  Alexandrite  Laser  Source 


DIAL  applications  require  precise  tuning  of  the  laser  transmitter's 
wavelength  output.  The  wavelength  needs  to  be  controlled  to  coincide  with 
narrow  absorption  lines  that  correspond  to  the  molecular  resonances  to 
be  studied.  Alexandrite  is  the  only  solid-state  laser  medium  that  is 
continuously  tunable  over  a  broad  spectrum  in  the  near-infrared,  720-  to 
790-nra  (Walling  et  al.,  1980). 

Allied  Chemical  Corporation  has  patented  a  process  to  grow  pure 
alexandrite  crystal  rods  synthetically.  This  material  exhibits  similar 
high  efficiency,  high  average  power,  and  short  pulse  duration  as  other 
solid-state  laser  media,  e.g.,  YAG  and  ruby.  In  essence,  the  design  of 
an  alexandrite  laser  resembles  that  used  for  YAG  laser  rods  (Koechner, 
1976).  Thus,  the  reliability,  ruggedness,  compactness,  and  ease  of 
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FIGURE  2  ONE-WAY  TRANSMISSION  SPECTRUM  FROM  13,600  TO  13,630  cm'1 
THROUGH  300  m  OF  1962  U.S.  STANDARD  ATMOSPHERE  AT  284  K 

All  absorption  lines  in  this  region  are  due  to  water  vapor. 
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13,600 


13.630 


FIGURE  4 


WAVENUMBER - cm-1 

ONE-WAY  TRANSMISSION  SPECTRUM  THROUGH  300  m 
OF  1962  U.S.  STANDARD  ATMOSPHERE  AT  314  K 


TRANSMISSION 


FIGURE  t  TRANSMISSION  LINE  PROFILES  AT  THREE  TEMPERATURES 
FOR  WATER  VAPOR  LINE  CENTERED  AT  13,602.63  cm'1 
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13.628.70 

FIGURE  6 


WAVENUMBER  - cm'1 

TRANSMISSION  LINE  PROFILES  AT  THREE  TEMPERATURES 
FOR  WATER  VAPOR  LINE  CENTERED  AT  13,628  88  cm'1 


13,629.00 
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254 

FIGURE  7 


TEMPERATURE  - deg  K 

TEMPERATURE  DEPENDENCE  OF  TRANSMISSION  AT  LINE  CENTER 
FOR  THREE  WAVENUMBERS  USED  WITH  THE  DIAL  WATER  VAPOR 
AND  TEMPERATURE  TECHNIQUE 
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operaci'n  (excellent  thermal  conductivity  allows  room  temperature  opera¬ 
tion;  ot  a  solid-state  laser  is  available  with  a  spectral  output  that 
is  generallv  accessible  only  with  dye  lasers  (which  are  usually  pumped  bv 
solid-state  lasers).  The  alexandrite  laser  medium  operates  as  ruby  ex¬ 
cept  that  the  terminal  lasing  level  is  a  band  of  vibrational  energy 
states;  these  yield  its  primary  tuning  range  (Walling,  1982).  This 
tuning  range  can  be  continuously  and  rapidly  scanned.  Figure  8  shows 
the  typical  output  energy  available  as  a  function  of  the  laser's  main 
tuning  range. 


FIGURE  8  TYPICAL  ENERGY  OUTPUT  OF  THE  ALEXANDRITE  LASER 
(AT  ROOM  TEMPERATURE  OPERATION)  AS  A  FUNCTION 
OF  OUTPUT  WAVELENGTH  TUNING 


Figure  9  is  a  general  schematic  containing  the  major  element,  of  a 
pulsed  alexandrite  laser.  A  pockels  cell  performs  the  Q-switchins  opera¬ 
tion  necessary  for  narrow  pulsed  outputs.  The  birefringent  tuner  controls 
the  coarse  continuous  wavelength  tuning.  Intraca.  '  y  etj'n  f.n  regulate 
fine  tuning  of  the  wavelength  and  the  linewidth  o:  the  output  pulse 
(Sam  and  Roukard,  1982). 

Table  4  lists  the  specifications  of  an  alexandrite  laser  that  are 
commercially  available.  These  capabilities  justify  the  use  of  the 
alexandrite  laser  as  a  DIAL  transmitter.  For  most  of  the  items,  improve¬ 
ment  through  an  arrangement  with  the  manufacturer  during  construction 
of  the  laser  can  be  achieved. 
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FIGURE  9  GENERAL  SCHEMATIC  OF  A  PULSED,  WAVELENGTH 
TUNABLE  ALEXANDRITE  LASER 


Table  4 

ALEXANDRITE  LASER  OUTPUT  SPECIFICATIONS 


Wavelength  (ran)  720-790 

Energy  (mj /pulse)  100 

Maximum  repetition  rate  (Hz)  10 

Linewidth  (ran)  .005 

Beam  divergence  (mrad)  2 
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SUMMARY  AND  RECOMMENDATIONS 


A.  Results 


A  study  and  review  of  experimental  systems  using  the  DIAL  and  Raman 
remote  systems  techniques  indicate  that  a  DIAL  system  should  be  constructed 
to  obtain  ground-based  vertical  water  vapor  and  temperature  measurements 
for  the  removal  of  certain  meteorological  system  errors  from  astronomical 
observations.  This  DIAL  technique  relies  on  a  laser  system  that  trans¬ 
mits  high  energy  pulses  at  three  specific  wavelengths  in  the  infrared 
spectrum. 

Such  DIAL  systems  have  typically  used  dye  lasers  (optically  pumped 
from  another  laser  source,  e.g.,  YAG)  for  study  in  the  700-  to  800-nm 
wavelength  region.  However,  the  alexandrite  solid-state  material,  which 
lases  with  continuous  tunability  in  the  same  spectral  region,  offers 
several  advantages: 

•  Chemical  stability  and  physical  durability. 

•  Simple  Q-switching  for  narrow  pulse  widths  because  of  strong 
polarized  output  from  the  crystal. 

•  Improved  output  energies  and  efficiency. 

•  Tuning  by  a  low-loss  quartz  birefringent  crystal  and  intra¬ 
cavity  etalons. 

Mass  production  of  alexandrite  rods  has  only  recently  been  accom¬ 
plished,  and  they  are  available  only  from  a  single  supplier.  (An 
alexandrite  laser  system  can  be  purchased  from  the  Military  Laser  Systems 
Division  of  Allied  Chemical  Corporation.)  Experiments  performed  over 
the  last  two  years  at  a  number  of  research  institutions  indicate  that  it 
will  be  a  significant  laser  source. 

The  alexandrite  laser  can  be  built  into  an  optical  and  data  collec¬ 
tion  system  that  will  collect  differential  absorption  data  and  store 
the  results  on  magnetic  tape  for  later  analysis.  Temperature  and  humidity 
values  should  also  be  collected  simultaneously  by  the  data  system  from 
in-situ  sensors  near  the  viewing  path.  Processing  of  the  correction 
algorithms  from  atmospheric  data  does  not  require  a  "real-time"  capa¬ 
bility. 

The  absorption  data  obtained  from  the  DIAL  system  will  be  trans¬ 
lated  into  estimates  of  humidity  and  temperature.  Correlative  compari¬ 
sons  and  checking  of  the  validity  of  results  can  then  be  made  with  the 
in-situ  sensors. 
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The  DIAL  technique  has  been  successfully  applied  in  the  724-nm  wave¬ 
length  region  to  obtain  absolute  humidity  measurements.  Also,  DIAL 
remote  temperature  measurement  systems  with  infrared  laser  sources  have 
used  oxygen  lines  near  760  run  (Mason,  1975;  Schwemmer  and  Wilkerson, 

1979).  However,  to  provide  simultaneous  humidity  and  temperature  mea¬ 
surements,  the  temperature  dependence  of  water  vapor  lines  near  734  nm 
needs  to  be  utilized.  The  suggested  DIAL  temperature  measurement  tech¬ 
nique  uses  the  line  strength  dependence  of  the  absorption  lines  centered 
near  735  nm  and  733  nm.  The  absorption  line  at  735  nm  deepens  with  in¬ 
creasing  temperature;  the  line  at  733  nm  absorbs  less  with  increasing 
temperature.  By  measuring  this  differential  change  in  absorption,  an 
estimate  of  the  absolute  temperature  can  be  calculated.  Thus,  to  obtain 
simultaneous  temperature  and  humidity  measurements,  a  three-wavelength 
DIAL  system  is  needed;  a  two-wavelength  system  can  provide  only  tempera¬ 
ture  or  humidity  results. 

Figure  10  shows  the  optical  components  needed  for  the  DIAL  system. 

As  shown,  three  laser  sources  will  fire  sequentially  with  200  ys  separa¬ 
tion  between  pulses.  The  short  inter-pulse  time  minimizes  the  temporal 
difference  in  the  atmosphere  that  each  pulse  could  encounter.  The  three 
beams  are  combined  with  beam  splitters  with  a  small  portion  of  the  energy 
removed  and  measured  by  a  photo-multiplier  tube  (PMT)  as  a  transmit  energy 
monitor.  The  combined  beams  are  then  emitted  through  the  center  of  the 
receiver  telescope,  allowing  common  optical  alignment  between  the  trans¬ 
mitted  and  received  beams.  The  receiver's  field  of  view  is  slightly 
larger  than  the  transmitted  beam  divergence  so  that  a  slight  error  in 
alignment  can  be  allowed.  Direction  of  the  optical  beams  is  controlled 
vertically  by  a  rotating  mirror  and  horizontally  by  a  dome  that  can  be 
rotated. 

The  backscattered  signals  are  collected  and  focused  onto  a  second 
PMT  that  measures  the  received  signal  strength.  Signal  amplification  in 
PMTs  is  accomplished  by  the  conversion  of  photons  to  electrons  by  multiple 
grids  within  the  unit.  An  optical  conversion  efficiency  of  20  percent 
can  easily  be  achieved.  Typically,  the  electronic  noise  level  output 
of  the  detector  is  lower  than  the  background  optical  signal  due  to  solar 
scattering  during  daytime  operation.  (The  use  of  narrow  bandwidth  optical 
interference  filters  centered  at  the  transmitted  wavelengths  minimize 
this  background  radiation  level.) 

The  received  signal  is  sampled  and  digitized  [with  a  high-speed 
(20  MHz)  10-bit  analog-to-digital  (A-to-D)  converter]  and  stored  in  the 
local  data  collection  system.  The  receiver  electronics,  consisting  of 
a  linear  amplifier  and  A-to-D  converter,  produces  a  representation  of 
the  returns  with  an  accuracy  of  .1  percent  over  a  range  of  three  decades 
in  magnitude.  A  sampling  interval  of  50  nsec  allows  a  range  resolution 
of  7.5  m  along  the  optical  path. 

By  measuring  the  time  delay  from  the  laser  firing  until  the  pulse 
reaches  a  region  of  interest  and  returns,  the  specific  backscattered 
strength  of  one  of  the  transmitter  wavelengths  can  be  obtained.  Even  though 
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a  single  receiver  channel  is  used  to  measure  data  for  all  three  wave¬ 
lengths,  specific  wavelength  backscatter  in  the  received  waveform  can 
be  selected  by  the  controlled  inter-pulse  time  of  the  laser  transmitters. 
With  this  information  the  DIAL  processing  equations  for  temperature  and 
water  vapor  density  can  be  utilized  during  off-line  processing,  and  a 
range  profile  can  be  estimated. 


B.  Recommendations 


Two  major  atmospheric  measurement  problems  for  compensation  of  dis¬ 
tortion  in  ground-based  astronomical  imaging  are  the  vertical  profile  of 
water  vapor  and  the  distribution  of  temperature  from  ground  level  to  a 
height  of  a  few  kilometers.  The  suggested  remote  sensing  lidar  system 
is  designed  specifically  to  obtain  simultaneous  data  of  these  two  quanti¬ 
ties.  It  consists  of  a  three-wavelength  laser  source  and  subsequent 
DIAL  analysis  of  the  range-resolved  backscattered  signals. 

The  technique  for  temperature  measurement  employs  all  three  wave¬ 
lengths,  (  two  of  which  are  tuned  to  water  vapor  absorption  lines  (near 
734  nm)  and  one  nonabsorbing  reference  wavelength.  Water  vapor  ontent 
can  be  obtained  with  a  DIAL  analysis  of  either  of  the  two  absorbing 
wavelengths  and  the  reference  wavelength. 

Alexandrite  should  be  employed  as  the  lasing  medium  for  the  lidar 
system,  allowing  ease  of  operation,  solid-state  ruggedness,  and  reli¬ 
ability.  The  laser  is  continuously  tunable  from  720  nm  to  790  nm, 
covering  an  important  spectral  region  for  water  vapor  absorption.  (For 
the  three-wavelength  metl  jd  at  734  nm,  two  of  the  wavelengths  should  be 
controlled  to  within  .005  nm  and  the  other  to  within  .1  nm. )  This  sys¬ 
tem  also  has  the  capability  to  study  other  techniques  in  the  infrared, 
including  water  vapor  absorption  measurements  at  724  nm  and  temperature 
and  density  measurements  at  760  nm  using  oxygen  absorption  lines. 

Before  a  system  can  be  directly  applied  to  studying  astronomical 
refraction  at  the  Naval  Observatory  in  Washington,  D.C.,  a  breadboard 
should  be  constructed  with  two  alexandrite  lasers.  The  primary  goals 
for  the  breadboard  system  would  be:  (1)  characterization  and  operational 
experience  with  the  alexandrite  lasers,  and  (2)  application  and  subse¬ 
quent  analysis  of  the  DIAL  technique  for  obtaining  water  vapor  profiles. 

A  third  laser  and  conical-scanning  optics  would  then  be  added  to  include 
simultaneous  temperature  distribution  and  isothermal  tilt  measurements. 
The  results  of  the  complete  breadboard  project  would  demonstrate  the 
functional  aspects  of  the  alexandrite  DIAL  system  for  the  remote  sensing 
of  water  vapor  and  temperature. 
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Abstract:  The  need  for  a  real  time  (or  near  real  time)  detailed  know¬ 
ledge  of  the  atmospheric  structure  necessary  for  determining  astronom¬ 
ical  refraction  in  general  is  exhibited.  Examples  involving  water 
vapor  and  isopycnic  tilts  are  given. 

The  possibility  of  determining  the  necessary  parameters  by  means  of 
active  atmospheric  probing  using  Light  Detection  and  Ranging  (Lidar) 
methods  is  reviewed.  A  brief  discussion  of  Raman  versus  Differential 
Absorption  Lidar  (DIAL)  techniques  is  included.  Recent  advances  in 
applicable  lasing  materials  such  as  the  Alexandrite  crystal  offer 
interesting  possibilities. 

At  the  present  time  a  conceptual  design  study  is  underway.  This  is 
intended  to  culminate  in  the  construction  of  a  system  in  1982-83. 

INTRODUCTION 

The  most  exasperating  characteristic  of  astronomical  refraction 
is  the  fact  that  while  it  is  a  relatively  straight-forward  matter  to 
account  for  nearly  99.9 %  of  the  effect,  it  is  currently  essentially 
impossible  to  allow  for  the  remaining  0.1%  with  any  assurance  of  suc¬ 
cess.  We  have  however,  reached  the  point  where  that  remaining  0.1%  is 
important,  and  especially  so  if  any  part  of  it  is  systematic.  This 
situation  is  not  caused  by  a  lack  of  excellent  theories  of  refraction. 

We  have  the  Pulkovo  Tables  (1956)  and  Garfinkel's  (1967)  polytropic 
theory  to  name  but  two  of  the  prominent  methods  available.  The  problem 
is  that  it  is  just  too  much  to  expect  any  theory,  even  one  with  many 
adjustable  parameters,  to  account  for  the  minute  by  minute  or  even 
hourly  variations  in  the  dynamic  planetary  boundary  layer.  On  the  whole 
it  is  noteworthy  how  well  the  theories  succeed  since  they  are  essentially 
static  representations  of  the  atmosphere. 

WATER  VAPOR 

Consider  the  effect  of  water  vapor.  It  changes  both  the  refractivity 
and  the  gas  constant  so  as  to  decrease  the  refraction.  Difficulties 
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arise  when  one  uses  only  a  ground  level  value  of  the  humidity  without 
any  knowledge  of  the  actual  three  dimensional  distribution  of  the  water 
vapor.  Merely  using  the  observed  ground  level  value  of  the  ratio  of 
the  partial  pressure  of  water  vapor  to  the  total  pressure,  Pw/Pt,  is 
equivalent  to  assuming  that  this  ratio  is  a  constant  throughout  the 
atmosphere.  Since  humidity  falls  rapidly  with  height  this  procedure 
would  generally  overestimate  the  effect  of  the  water  vapor. 

In  order  to  illustrate  this  type  of  problem  the  changes  in  re¬ 
fraction  (with  respect  to  dry  air)  caused  by  various  partial  pressures 
of  water  were  computed  using  the  Pulkovo  Tables  (IV  ED.),  and  also  by 
means  of  numerical  integrations  of  various  water  vapor  profiles.  The 
differences  in  the  changes  using  these  two  methods  were  then  computed 
and  tabulated.  For  a  total  ground  level  pressure  near  760  mm  Hg, 
with  10  mm  being  due  to  water  vapor,  and  for  temperatures  in  the  range 
of  273°  to  293°  K,  it  was  found  that  good  agreement  occurred  when  the 
water  vapor  was  assumed  to  decline  linearly  with  height,  reaching  zero 
at  seven  kilometers.  Table  I  illustrates  these  differences.  Change 
(Integration)  -  Change  (Pulkovo),  for  that  profile.  Evidently  if  one 
happened  to  be  observing  under  just  the  right  conditions,  a  near  perfect 
allowance  for  water  vapor  would  be  made.  However,  as  Table  I  shows,  the 
differences  become  greater  at  lower  total  pressures  or  higher  partial 
pressures  of  water  vapor.  There  is,  in  addition,  no  guarantee  that  the 
true  profile  is  that  which  was  assumed  here. 


TABLE  I 

Change  (Integration)  -  Change  (Pulkovo) 


\T°K 

\  273 

283 

293 

Pt  mm 

Pw  -  10 

Pw  -  30 

Pw  -  10 

Pw  -  30 
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t— * 
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Pw  =>  30 

ff 

ft 

f* 

If 

»» 
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760 

-.001 

-.014 

-.001 

-.012 

-.002 

-.013 

640 

+.017 

+  .045 

+  .018 

+  .042 

+  .017 

+  .042 

560 

+.032 

+.084 

+  .030 

+  .081 

+  .029 

+  .079 

ISOPYCNIC  TILT 


A  second  contributor  to  the  types  of  error  considered  here  is  any 
systematic  tilt  of  the  assumed  isopycnics.  One  of  us,  Hughes  (1979), 
has  done  considerable  theoretical  work  on  the  generation  of  such  tilts 
by  urban  heat  islands.  Without  going  into  great  detail  we  present  here 
two  cross  sections  through  the  density  perturbation  field  caused  by  a 
realistic  heat  island.  Cross  section  A  lies  in  the  plane  of  the  meridian 
of  a  hypothetical  observatory,  while  section  B  lies  in  the  plane  of  the 
corresponding  prime  vertical.  Table  II  shows  the  values  of  the  tilt 
angles  in  these  planes.  The  observatory  is  supposed  to  be  located  some 
6  km  southeast  of  the  center  of  the  heat  island.  The  wind  is  assumed  to 
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be  from  the  west.  Values  are  given  up  *-o  a  height  of  500  m  and  for 
horizontal  distances  of  400  m  and  800  m  (negative  to  the  south  for 
section  A  and  to  the  west  for  section  B) . 


TABLE  II 
Isopycnic  Tilts 


Section 

A 

Section  B 

Meters 

-800 

-400 

0 

+400 
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-800 
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0 
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These  results  are  strictly  theoretical,  and  although  the  order  of  magni¬ 
tude  of  the  tilts  is  reasonable,  there  is  no  guarantee  that  the  true 
tilt  is  that  which  is  shown  here. 

Thus  for  the  two  phenomena  considered  above  we  are  presently  forced 
into  what  are  essentially  ad  hoc  assumptions  and  procedures.  This  un¬ 
satisfactory  situation  can  only  be  improved  by  injecting  additional, 
real  time,  three  dimensional  atmospheric  data  into  the  various  algori¬ 
thms  used  for  treating  observational  material. 

LIDAR 

One  method  of  generating  the  required  data  is  with  a  laser  detection 
and  ranging  (Lidar)  system.  The  approach  is  similar  to  conventional 
radar,  except  the  interactions  occur  at  optical  frequencies.  With  the 
laser  transmitter  and  receiver  systems  positioned  at  a  common  location, 
a  telescopic  optical  receiver  measures  the  amount  of  backscatter  from 
a  short  transmitted  laser  pulse.  Ranging  is  performed  by  calculating 
the  time  for  the  pulse  to  travel  the  distance  to  a  region  of  interest 
plus  the  time  for  the  backscattered  pulse  to  return.  Thus,  by  sampling 
the  time-varying  received  optical  intensity  after  a  specific  delay, 
study  of  a  restricted  location  along  the  path  can  be  accomplished. 

By  using  models  that  describe  the  absorption  and  scattering  properties 
of  the  atmosphere  at  optical  frequencies,  information  about  major  and 
trace  constituents  can  be  obtained  through  "inverting"  the  measurements 
of  backscattered  intensity.  Depending  on  the  realm  of  study,  the  fre¬ 
quency  of  the  transmitter,  the  frequency  response  of  the  receiving  system, 
the  optics  design,  and  the  data  processing  can  be  selected  to  provide  a 
large  area  of  coverage  and  rapid  measurement  capabilities.  See  Hinkley 
(1976)  for  an  excellent  review  of  the  field. 
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DIAL  AND  RAMAN  LIDAR  SYSTEMS 

Two  different  Lidar  techniques  will  be  discussed  in  this  paper:  (1) 
Differential  Absorption  Lidar  (DIAL)  and  (2)  Raman  Lidar.  Each  of 
these  are  candidates  as  system  approaches  to  provide  vertical  profiles 
of  atmospheric  water  vapor  content  and  temperature.  However,  the 
scattering  and  absorption  basis  for  the  measurements  are  distinct  in 
each  system. 

The  laser  transmitter  of  the  DIAL  system  sends  pulses  at  two  wave¬ 
lengths.  One  wavelength  is  selected  to  coincide  with  an  absorbing  line 
of  the  molecular  species  of  interest.  The  other  near-by  wavelength  is 
chosen  for  negligible  absorption.  The  receiving  system  is  tuned  to 
measure  the  elastic  Mie  and  Rayleigh  backscatter  at  the  same  wave¬ 
lengths  as  were  transmitted.  By  forming  the  difference  of  the  returned 
intensities  at  the  absorbed  and  unabsorbed  wavelengths,  a  time-varying 
(and  thus,  range  resolved)  concentration  of  the  species  is  estimated. 

This  technique  has  also  been  called  Differential  Absorption  of  Scattered 
Energy  (DASE)  and  Differential  Absorption  and  Scattering  (DAS). 

The  transmitter  of  a  Raman  Lidar  system  sends  a  pulse  of  a  single 
wavelength.  The  receiver  system  is  tuned  to  the  wavelengths  that  are 
uniquely  Raman-shifted  (from  the  transmitted  wavelength)  by  the  mole¬ 
cular  species  to  be  studied.  By  using  the  Raman  cross-section  of  the 
species,  an  analysis  of  the  magnitude  of  the  received  intensity  can  be 
translated  into  a  concentration  profile. 

WATER  VAPOR  MEASUREMENT  SYSTEMS 

Both  Raman  and  DIAL  systems  are  used  successfully  to  obtain  water 
vapor  measurements.  Only  those  experiments  that  have  supplied  vertical 
profiles  are  listed  in  the  following  tables;  those  systems  that  supply 
a  column  integrated  measurement  or  average  along  a  horizontal  path  are 
not  included. 

Table  III  shows  the  range  of  Raman  systems  that  have  been  designed 
specifically  for  water  vapor  profiles.  All  of  these  systems  fundamentally 
measure  the  ratio  of  the  number  density  of  H2O  to  that  of  N£>  Thus, 
only  a  water  vapor  mixing  ratio  is  estimated.  This  can  be  translated 
into  a  water  vapor  profile  by  assuming  a  density  profile  for  N2  with 
accompanying  error. 

Table  IV  lists  recent  water  vapor  profiling  DIAL  systems.  Less 
powerful  laser  transmitters  are  typically  needed  than  those  used  in 
Raman-based  systems  because  the  Rayleigh-Mie  backscatter  cross-sections 
are  many  orders  of  magnitude  larger  than  Raman  cross-sections. 
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TABLE  III 

Raman  Lidar  Systems  for  Vertical  Water  Vapor  Profiling 


Cooney  (1971) 

Transmitter:  doubled  Ruby  at  347.15nm.  Receiver:  397.5  nm  from 
H2O  and  377.7  nm  from  Approach:  water  vapor  profile  obtained 

by  normalizing  H2O  Raman  return  with  N2  return,  and  assuming  an 
N2  atmospheric  content.  13  percent  accuracy.  Night  measurement. 
Range:  300  m  -  2600  m  with  100  m  intervals. 


Melfi  (1972) 

Transmitter:  doubled  Ruby  at  347.15  nm.  Receiver:  397.5  nm  from 
H2O  and  377.7  nm  from  N2.  Approach:  water  vapor  profile  from 
H2O  return  normalized  by  N2  return.  Night  measurement.  Range: 

200  m  -  2500  m  with  50  m  intervals. 


Pourny,  Renault,  and  Orszag  (1979) 

Transmitter:  doubled  ruby  at  347.15  nm.  Receiver:  397.6  nm  from 
H2O  and  377.7  nm  from  N2.  Approach:  water  vapor  mixing  ratio  fron 
H2O  return  normalized  by  N2  return.  Night  measurement.  15  percent : 
accuracy.  Range:  300  m  -  1800  m  with  30  m  intervals. 


Renault,  Pourny,  and  Capitini  (1980) 

Transmitter:  quadrupled  YAG  at  266  nm.  Receiver:  277.5  nm  from 
O2.  283.6  nm  from  N2,  and  294.6  nm  from  H2O.  Approach:  water 

vapor  mixing  ratio  from  the  H2O  return  is  normalized  by  the  N2 
return.  Estimation  of  ozone  concentration  is  used  to  compensate 
for  heavy  ozone  attenuation  of  the  signals  in  this  wavelength 
region.  10  percent  accuracy  at  500  ra  with  30  m  intervals.  Day¬ 
time  operation.  Range:  150  m  -  950  m  with  30  m  intervals. 
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TABLE  IV 

DIAL  Systems  for  Vertical  Water  Vapor  Profiling 


Browell,  Wilkerson,  and  Mcllrath  (1979) 

Transmitter:  ruby  at  694.3  nm  and  dye  (pumped  by  ruby)  at  724.3 
nm.  Receiver:  724.3  nm  (water  vapor  absorbing  wavelength)  and 
694.3  nm.  Approach:  difference  in  returned  signals  is  analyzed 
for  H2O  density  profile  (mol/cm^).  Equipment  includes  an  absorption 
cell  that  measures  the  cross-section  for  H2O  at  the  dye  laser  out¬ 
put  for  each  pulse.  9  percent  accuracy.  Night  operation.  Range: 
100  m  -  3000  m  with  100  -  180  m  intervals. 


Cahen,  Pelon,  Flamant,  Lefere,  Chanin,  and  Megie  (1980);  also,  Cahen 
and  Lesne  (1981) 

Transmitter:  Nd:  YAG  at  532  nm  and  dye  (pumped  by  Nd:YAG)  at  723 
run.  Receiver:  723  nm  (absorbed  by  water  vapor)  and  532  nm 
Approach:  difference  returns  to  obtain  water  vapor  profiles. 

Wavelength  and  line-width  are  servo-controlled.  10  percent  accuracy. 
Daytime  operation.  Range:  8000  m  with  30  m  intervals. 


Werner  and  Herrmann  (1981) 

Transmitter:  ruby  at  694.28  nm  and  temperature  varied  ruby  from 
694.2  nm  to  694.28  nm.  Receiver:  694.215  nm  (absorbed  by  water 
vapor)  and  o94.237  nm.  Approach:  difference  returns  to  obtain 
water  vap  r  density  profiles.  Accuracy  of  1  mbar  in  absolute 
humidity.  Daytime  operation.  Range:  450  m  -  1500  m  with  100  m 
intervals. 


TEMPERATURE  PROFILING  SYSTEMS 

Table  V  contains  a  list  of  Raman  temperature  profiling  systems. 
Common  to  these  experiments  is  the  technique  of  measuring  the  temperature 
dependence  in  the  magnitude  of  the  Raman  spectrum  of  N2  and  O2.  The 
lines  in  the  Raman  spectrum  have  differential  changes  in  magnitude  wi hh 
respect  to  temperature.  Thus,  by  comparing  the  change  in  two  bands, 
temperature  value  can  be  deduced. 


There  is  a  lack  of  DIAL  vertical  temperature  profiling  systems. 
However,  work  is  presently  being  performed  with  system  techniques  that 
show  promise.  One  proposal  is  to  use  "Mason's  Method"  (Mason,  1975- 
Schwemmer  and  Wilkerson,  1979)  with  the  absorption  spectra  of  O2  in  the 
region  of  760  nm.  This  wavelength  is  close  to  the  wavelengths  selected 
in  successful  DIAL  water  vapor  measurements  (720  nm) .  Thus  a  single 
system  could  supply  temperature  as  well  as  water  vapor  density. 
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TABLE  V 

Raman  Lidar  Systems  for  Vertical  Temperature  Profiling 


Cooney  and  Pina  (1976) 

Transmitter:  ruby  at  694.3  nm.  Receiver:  689.0  nm  and  691.2  nm 
Approach:  measure  the  difference  in  the  magnitude  of  the  Raman 
spectrum  of  O2  and  N2  in  tvo  bands  which  vary  with  temperature. 
Accuracy  is  5  degrees  C.  Nighttime  operation.  Range:  600  m  - 
1100  m  with  300  m  intervals. 


Gill,  Geller,  Farina,  and  Cooney  (1979) 

Transmitter:  ruby  at  694-. 3  nm.  Receiver:  691.6  nm  and  689.0  nm. 
Approach:  same  as  above.  Accuracy  is  .85  degrees  C.  Range: 

1300  m  -  2300  m  with  75  m  intervals. 


A  COMBINATION  SYSTEM 

The  combined  DIAL  technique  mentioned  above  is  now  being  researched. 
Typically,  dye  lasers  have  been  used  to  cover  the  650  nm  -  790  nm 
region  as  transmitters  for  such  simultaneous  temperature  and  humidity 
measurement  system  studies.  Besides  covering  the  important  absorption 
regions  for  H2O  and  02*  they  are  easily  tuned  to  specific  wavelengths. 
However,  solid-state  lasers  have  demonstrated  an  ease  of  operation  and 
reliability  that  is  important  for  day-to-day  meteorological  observations. 
Recently,  a  solid-state  lasing  material.  Alexandrite,  was  developed  for 
laser  output  from  700  nm  -  815  nm.  Thus,  the  characteristics  of  this 
material  are  being  included  within  new  DIAL  system  designs. 

CONCLUSION 

It  appears  that  Lidar  offers  real  possibilities  with  regard  to  pro¬ 
viding  input  data  for  use  in  refining  our  current  allowances  for  refrac¬ 
tion.  Water  vapor  profiling  presents  no  problem,  and  temperature  measure¬ 
ments  certainly  are  feasible.  Since  for  the  small  perturbations  consider¬ 
ed  here,  the  Boussinesq  approximation  is  valid,  it  follows  that  isopycnics 
follow  isotherms.  Thus  temperature  measurements  are  equivalent  to 
density  measurements,  and  with,  for  example,  a  conically  steered  laser 
beam,  the  isopycnic  tilts  can  be  determined. 

If  present  plans  hold,  a  system  to  accomplish  these  ends  will  be 
in  operation  within  two  years.  The  results  should  be  of  considerable 
interest. 
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